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ABSTRACT 
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O , 

The first results from a near-contemporaneous optical and infrared spectro- 
scopic observing program designed to probe the detailed density structure of 



(N 



classical Be circumstellar disks are presented. We report the discovery of asym- 



> 



metrical infrared emission lines of He I, I, Fe II, and the Brackett, Paschen, 
and Pfund series lines of H I which exhibit an opposite V/R orientation (V > R) 
to that observed for the optical Balmer Ha line (V < R) in the classical Be star 
( Tau. We interpret these data as evidence that the density wave which charac- 
terizes ( Tau's disk has a significantly different average azimuthal morphology in 
the inner disk region as compared to the outer disk region. A follow-up multi- 
wavelength observational campaign to trace the temporal evolution of these line 



O 

Pl,: 

6 

vidual (C Tau) — stars: individual (48 Librae) 
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profile morphologies, along with detailed theoretical modeling, is suggested to 
test this hypothesis. 
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Introduction 



Both observational and theoretical evidence indicates that classical Be stars are near 



main sequence B-type stars with gaseous circumstellar disks (e.g. see review by lPorter fe Rivinius 
20031 ); however, the mechanism(s) responsible for creating these disks are not well under- 
stood. Optical spectroscopic studies have uncover ed a myriad of differe nt line profile mor- 
phologies commonly observed in class ical Be stars ( Hanuschik et aD 1996 ), i ncluding double 



peake d asymmetrical emission lines ([McLaughlin! Il96ll ; iDachs et al.l 



1992; Hanuschik et al. 



19951 ). The ratio of the violet (V) to red (R) flux within individual lines, hereafter referred 
to as a V/R ratio, is observed to vary on q uasi-periodic ti mescales of several to ten years, 
as documented for a large set of Be stars by lOkazakil (119971 ). While numer ous theories have 



been p roposed to explain this observed phenomenon, as summarized within lHanuschik et al. 
( 119951 ) . these features are most com monly interpre t ed as evidence of one-armed densit y 



waves within these circumstellar disks (10kazakilll99ll ; lHanuschik et al.lll995t IOkazakilll9971 ). 
No te however that t he pr esence of phase lags between polarimetric and Ha observations has 
led iMcDavid et al.l (120001 ) to sugge st these waves might be more spiral -like than previously 
thought. Interferometric studies (jVakili et a.l 1 19981 ; iBerio et al.l 119991 ) have also provided 
evidence that Be disks deviate from axisymmetry, namely they have provided evidence of 
prograde one-armed oscillations in several disks. 

It is not uncommon to observe multiple emission lines in a spectrum of a classical Be 
star that exhibit a V/R asymmetr y, and with f ew exceptions the same orientation of this 



19821) and be tween different lines in the hydrogen Balmer series (IKogure fe Suzukil 



and Fe II lines f 


Sle 


ttcbak 


(Kogure & Suzuki 


1984 





Baadd (119851 ) also suggest ed that a phase lag exis t s bet ween the Ha and hydrogen Pa 15-18 
lines of 7 Cas presented by lChalabaev fe Maillardl (119831 ); however, we believe this suggestion 
is rather dubious given the poor resolution of these data and clear evidence of contamination 
from nearby lines. 

Few moderate resol ution infrared ( I R) sp ectro scopic studies of class ical Be stars have 
been performed to date. IClark fe Steele! (120001 ) and lSteele fe Clark! (120011 ) have investigated 
a large set of Be stars in the H and K band at sufficient resolution to begin to search for 
detailed line profile effects. The first tenta tive evi d ence of IR V/R reversal s, analogous to 
the op tical V/R phase lags summarized by iBaadei (119851 ) . was presented in IClark fe Steele 
( 120001 ). who noted that the Brackett 7 and He I 2.058/iin lines had opposite orientations in 
one of their targets. 



In this Letter, we present the initial results from a near-contemporaneous, optical and 
IR spectroscopic observing campaign designed to provide insight into the detailed density 
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structure of classical Be circumstellar disks, hence diagnose the mechanisms responsible for 
creating these disks. Detailed modeling of these data, as well as follow-up observations, will 
be presented in subsequent papers. In Section 2 we outline the observational data to be 
discussed in this work. An analysis and discussion of these data are presented in Section 3. 
Finally, a summary of our main results and suggested follow-up work is presented in Section 
4. 



2. Observations 



The infrared (IR) spectroscopic data presented in this paper were obtained at the 3m 
NASA Infrared Tele scope Facility (IRTF), using the SpeX medium resolution spectrograph 
(IRayner et al.ll2003l ). We used a 0.3 x 15 arc-second slit to provide wavelength coverage from 
0.8-2.4/xm at R ~2000. Observations of AOV stars were obtained at ~15 minute intervals 
throughout our run, at a similar hour angle and airmass to our science targets, to serve 
as telluric standards. Further details regarding the tellu ric correction method commonly 
applied to Spe X data are presented in IVacca et al.l (120031 ). All observations were reduced 
using Spextool (jCushing et al.l I2004T ) . These data were taken during a period of extremely 
poor weather conditions on Mauna Kea, thus we have not attempted to perform a rigorous 
flux calibration. While these observing condit ions did introduce small artifacts in wavelength 
regimes outside of those presented here (see IWisniewskil 120051 ). we are fully confident that 
the line profile morphologies and line-to-continuum ratios of the data presented in this paper 
are not influenced by the presence of any residual telluric features. 

The optical spectroscopic data presented here were obtained at the lm telescope of the 
Ritter Observatory, using its fiber-fed echelle spectrograph. Nine non-overlapping orders, 
each of width 70 A, were obtained over the wavelength range 5300-6600A at R ~26000. 
Standard IRAF0 techniques were used to reduce these data, including procedures to apply 
bias, flat field, and wavelength calibration corr ections to all da t a. Fu rther details regarding 
the reduction of Ritter data can be found in iMorrison et al.l (119971 ). A summary of the 
optical and IR data to be discussed in this paper are presented in Table [TJ Note that we 
have not corrected any of our optical or IR spectral line data for underlying photospheric 
absorption components. 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Asso- 
ciation of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation. 
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Results and Discussion 



As demonstrated by our near-contemporaneous optical and IR observations of 48 Librae 
(see Figure [1]), double-peaked asymmetrical profiles in Be sp ectra often exhibit the same 
V/R orientation for all observed transitions. ISlettebakl (119821 ) reported phase lags between 
optical Fe II and hydrogen lines in several of his observations of Be stars; we found evidence 
of a similar phen omenon in our IR observation of the suggested classical Be star NGC 



2439:WBBel (see IWisniewskill2005t IWisniewski fe Biorkmanl 120061 ). 



We obtained two Ha spectra of the well studied classical Be star ( Tau (see Figure 
[2]) one month before and one month after our IR observations, and note that both Ha 
profiles exhibit triple-peaked profiles characterized by an overall V < R asymmetry. The 
V/R period is o ften several ye ars, so the line profiles are generally stable on timescales of 
several months (jokazaki 1997 ). We are confident that ( Tau exhibited a Ha profile similar 
to that shown in Figure [2] at the time of our IR obser yation, and was n ot in the midst of a 
cycle of irregular fluctuations of the type described by iGuo et al.l (119951 ). based on our long 
baseline of previous observations of the star. At the resolution of SpeX, ( Tau's IR spectrum 
shows many emission-lines characterized by a range of morphologies, including single-peaked 
profiles, asymmetrical double-peaked profiles, and asymmetrical double-peaked shell profiles. 
As shown in Figure [2J hydrogen Brackett transitions in these data exhibit shell-like profiles, 
which are asymmetrical with V > R. Other IR emission-lines which have asymmetrical 
profiles with a similar V/R orientation to the H I Brackett series include O I 0.8446, 1.129, 
1.317 /im; Fe II 1.086, 1.113 /im; He 1 1.083 /im; H I Paschen 7 and /3; and the H I Pfund lines 
at 2.382, 2.392, 2.403 /im. These data exhibit a fundamentally different V/R morphology 
than our near-contemporaneous Ha profile. 

As ( Tau is known to ex perience cyclic al V/R variations, likely due to the presence 
of a density wave in its disk (jOkazakil 119971 ). we suggest that the observed opposite V/R 
orientation of the IR lines and Ha can be explained by the presence of a density wave in the 
disk which has a significantly different average azimuthal morphology (e.g. a spiral wave), 
combined with different loci for line formation. This idea is illustrated in Figure [3J where we 
show schematically the disk of ( Tau with a one-armed spiral density perturbation: note that 
lighter shading denotes regions of density enhancement while darker shading denotes density 
decrements. The H I IR lines are much more optically thin than the optical lines, hence they 
form within the inner disk, while the optical H I lines originate from a much larger area. As 
a result of the different azimuthal properties of the disk with respect to the distance from the 
star, in certain phases of the V/R cycle the net density enhancement of the inner disk may 
occur in an opposite side (with respect to the line of sight) to the net density enhancement 
of the entire disk, thus producing line profiles with opposite V/R ratios. For instance, in 
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Figure [3] the V > R morphology of Br-12 is produced by a net density enhancement in the 
inner disk region which is moving towards the observer, and the V < R morphology of Ha is 
produced by a net density enhancement in the outer disk region which is moving away from 
the observer. 



Carciofi fc Bjorkmanl (120061 ) , ICarciofi et al.l (120061 ) , and ICarciofi et al.l (120071 ) outline a 
new Monte Carlo code for solving the radiative transfer in the circumstellar disks of Be 
stars. This fully three-dimensional non-LTE Monte Carlo code produces a self-consistent 
determination of the hydrogen level populations, electron temperature, and gas density of 
a circumstellar disk, given a prescription of the central star and assuming a static, nearly- 
Keplerian dis k with a given mass lo ss rate, truncation radius, and a prescription for the gas 



viscosity. In ICarciofi et a 



(200 5|), t h e aut hors used this code to study the circumstellar 



disk of C Tau. Following I Wood et al.l (119971 ) . they adopted a luminosity class of IV and a 
spectral type of B2.9 ± 0.4, which corresponds to an effective temperature of 19, 000K and 
a stellar radius of 5.6R & . Assuming an outer disk radius of 100 R± and a mass loss rate 
of 2.5 x 10" 11 M /yr the authors were able to fit the observed SED from visible to mid 
IR wavelengths, and optical (AA 3200-10500 A) spectropolarimetry. Using this model, we 
calculated intensity maps of Ha and Br-12 for ( Tau which, as seen in Figure HJ illustrate 
that Br-12 is formed in the inner disk (R < 15i2*) whereas Ha is produced over a much 
more extended region, thereby confirming our previous assertion that our optical versus 
IR H I data probe distinctly different disk regions. Future iterations of the model will 
incorporate complex density profiles, such as one-armed waves, which will enable us to test 
our exploratory discussion in a more quantitatve manner. 

We note that V/R reversals have been predic ted by at least one other theory in the 
refereed literature. IWaters fc Marlborough I (119921 ) predicted instances in which a ratio of 



R> V in a Ha line would reverse to a ratio of R< V in IR H I lines, due to line self-absorption 
in the inner region of edge-on disks. However, the presence of this reversal depends criti- 
cally on the existence of a large expansion component in the disk velocity at the location 
where the lines form, an unlikely sc enario given the very small radial velocities now believed 



to ch aracterize inner disk regions (jRivinius et al 
2003h . 



1999 



Hanuschik 2000; Porter Rivinius 



4. Summary and Future Work 

We have presented the first results of near-contemporaneous optical and IR spectroscopic 
observations of classical Be stars. We report the discovery of asymmetrical infrared emission 
lines (He I, O I, Fe II, H I Paschen, Brackett, and Pfund series) in ( Tau which have oppositely 
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oriented V/R ratios to that seen at Ha. We suggest these data indicate that the density wave 
which characterizes ( Tau's disk has a significantly different average azimuthal morphology in 
the inner disk region as compared to its outer disk region. This would also indicate that the 
Fe II lines are formed in the inner disk, in the same spatial location as the IR lines are formed. 
These features should tr avel through th e disk on the cyclical timescale generally reported for 
density waves in ( Tau (jOkazakilll997l ); hence, we would predict continued anti-correlation 
amongst these line profiles over the course of several observing seasons. We have begun 
a long-term optical and IR spectroscopic monitoring program to trace this evolution. As 
demonstrated by the initial results presented in this Letter, such an observational approach 
has the potential to provide powerful constraints to theoretical efforts to map the detailed 
density distribution of Be circumstellar disks, hence better understand the mechanism(s) 
which create these disks. 
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Fig. I. — Left: The Ha spectrum of 48 Librae observed on 2004 March 11 exhibits an asym- 
metrical profile characterized by V > R. Right: A near contemporaneous infrared spectrum 
obtained on 2004 March 6 also exhibits hydrogen Brackett lines having asymmetrical profiles 
with V > R. The agreement between the orientation of these line profile morphologies is the 
standard scenario expected in observations of classical Be stars. 
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Fig. 2. — Left: The Ha spectrum of ( Tau obtained on 2004 February 28 (left panel) clearly 
exhibits an asymmetrical triple-peaked profile characterized by V < R. Although not shown 
here, an optical spectrum obtained on 2004 April 12 exhibits the same morphology. Right: 
This IR spectrum of £ Tau obtained on 2004 March 5 is characterized by asymmetrical 
hydrogen Brackett lines having V > R. Although not picture here, additional He I, O I, Fe 
II, and H I (Paschen and Pfund series) emission lines in this spectrum exhibit this same V 
> R morphology. 
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Fig. 3. — A one-armed density perturbation, characterized by regions of density enhancement 
(lighter regions) and regions of density decrements (darker regions), is believed to reside in 
C Tau's disk based on observed V/R line profile variations. Our IR (left panel) and optical 
(right panel) spectroscopic data probe different physical locations within the ( Tau disk, 
hence probe different regions of the density perturbation located within the disk. 
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Fig. 4. — ICarciofi fc Bjorkmanl (120061 ) model images of the ( Tau disk at Ha (upper panel) 



and H I Br-12 (lower panel). The scale of these images is 40 stellar radii (~224 solar radii). 
These images clearly demonstrate that Ha is produced over a much larger disk area, whereas 
H I Br-12 is produced from a more inner disk region. 



Table 1: Summary of Observations 



Object 


Observatory 


Date 


Exposure Times 


( Tau 


I 


2004 March 5 


60 s. 


CTau 


R 


2004 February 28 


1200 s. 


CTau 


R 


2004 April 12 


1200 s. 


48 Lib 


I 


2004 March 6 


150 s. 


48 Lib 


R 


2004 March 11 


3600 s. 



Note. — Some properties of our observational data are summarized. The observatory abbreviation "I" 
corresponds to IRTF data and "R" corresponds to Ritter data. For our infrared observations, the listed 
exposure times are the total effective integration times of the science targets, and do not represent the 
individual exposure times at each of the telescope nod positions described in the text. 



